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We review a certain class of ("nearly") exactly solvable models of electronic spectrum of two-dimensional 
systems with fluctuations of short range order of "dielectric" (e.g. antiferromagnetic) or "superconducting" type, 
leading to the formation of anisotropic pseudogap state on certain parts of the Fermi surface. The models are 
based on recurrence procedure for one- and two-electron Green's functions which takes into account of all Feynman 
diagrams in perturbation series with the use of the approximate Ansatz for higher-order terms in this series. These 
models can be applied to calculation of spectral density, density of states and conductivity in the normal state, 
as well as to calculation of some properties of superconducting state. 



The model of "nearly - antiferromagnetic" 
Fermi-liquid is based upon the picture of well de- 
veloped fluctuations of AFM short range order 
in a wide region of the phase diagram. In this 
model the effective interaction of electrons with 
spin fluctuations is described via dynamic spin 
susceptibility x q (w), which is determined mainly 
from the fit to NMR experiments : 



Veff(<l,bj) =g 2 X q (w) 
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where g is coupling constant, ^-correlation length 
of spin fluctuations, Q = (ir/a, 7r/a)-vector of 
antiferromagnetic ordering in insulating phase, 
a; .^-characteristic frequency of spin fluctuations, 
a-lattice spacing. As dynamic spin susceptibil- 
ity x q (w) has peaks at the wave vectors around 
(tt/q, 7r/a) there appear "two types" of quasipar- 
ticles — "hot quasiparticles" with momenta in the 
vicinity of "hot spots" on the Fermi surface and 
"cold" quasiparticles with momenta on the other 
parts of the Fermi surface, e.g. around diagonals 
of the Brillouin zone \p x \ = \Py\ [0 • 

In the following we shall consider the case of 
high enough temperatures when irT oj s f which 
corresponds to the region of "weak pseudogap" 
0. In this case spin dynamics is irrelevant and 
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Figure 1. "Hot spots" model. Electronic states 
around the intersection points of the Fermi sur- 
face with magnetic Brillouin zone (shown by 
dashed lines) are strongly interacting with fluc- 
tuations of AFM short range order. 



we can limit ourselves to static approximation. 

We can greatly simplify all calculations if in- 
stead of ([!]) we use another form of model inter- 
action: 



VefM) = A 5 
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where A is some phenomenological parameter, 
determining the effective width of the pseudogap. 
In fact (0) is qualitatively similar to the static 
limit of ([j]) and differs from it very slightly in 
most important region of |q — Q| < 

The spectrum of "bare" (free) quasiparticles 
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Figure 2. Model of the Fermi surface for HTSC- 
cuprates. "Hot patches" are shown by thick lines 
of the width of 



can be taken in the form ffl : 

£ p = — 2t(cos p x a+ cos p y a)— 4i cos p x a cos p y a(3) 

where ^-nearest neighbor transfer integral, t - 
second nearest neighbor transfer integral on the 
square lattice, \i is chemical potential. 

Much simpler model of electronic spectrum as- 
sumes that the Fermi surface of two - dimen- 
sional system to have nesting ("hot") patches of 
finite angular size a in (0, it) and symmetric di- 
rections in the Brillouin zone, as shown in Fig. 2 
Q . Similar Fermi surface was observed in a num- 
ber of ARPES experiments on cuprate supercon- 
ductors. Here we assume that fluctuations inter- 
act only with electrons from the "hot" (nesting) 
patches of the Fermi surface, and scattering vec- 
tor is either Q x = ±2pp, Q y ~ or Q y = ±2pp, 
Qx — for incommensurate fluctuations, while 
Q = (7r/a,7r/a) for commensurate case. It is eas- 
ily seen that this scattering is in fact of one - 
dimensional nature. In this case non-trivial con- 
tributions of interaction with fluctuations appear 
only for electrons from "hot" patches, while elec- 
trons on "cold" parts of the Fermi surface remain 
free. 

These models can be solved exactly in the limit 
of infinite correlation length £ — > oo, using meth- 
ods developed in Refs. For the case of finite 
£ we can use an approximate Ansatz, proposed 
for one - dimensional case in Ref. and further 



developed for two - dimensional system in Refs. 
0,0. According to this Ansatz the contribution 
of an arbitrary diagram for electron self-energy of 
iV-th order in interaction (M) has the form: 
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where for the "hot spots" model £ 3 - = £ p +q and 

v j = K+qI + K+qI for odd and & = £ P and 
v j = \ v p\ + l^pT^ 01 even 3 ~ appropiate com- 
binations of velocity projections, determined by 
the "bare" spectrum (||). For the "hot patches" 
model = (— 1) J £ P , Vj = vp. Here rij is 
the number of interaction lines, enveloping j-th 
Green's function in a given diagram. In this case 
any diagram with intersecting interaction lines is 
actually equal to some diagram of the same or- 
der with noncrossing interaction lines. Thus in 
fact we can consider only diagrams with nonin- 
tersecting interaction lines, taking into account 
diagrams with intersecting lines introducing addi- 
tional combinatorial factors into interaction ver- 
tices. This method was used for one-dimensional 
model of the pseudogap state in Refs. 

As a result we obtain the following recursion 
relation for one-electron Green's function (con- 
tinuous fraction representation) 0,^,0: 



G 1 (s n £,p) — G 1 (e n ^ p ) — S 1 (e n ^ p ) 



£ fc (e n £ p ) = A 2 



v(k) 



Combinatorial factor: 

v(k) = k 



(5) 
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corresponds to the case of commensurate fluctu- 
ations with Q = (n /a, n /a) [B. For incommensu- 
rate case [El: 



v(k) 
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for odd k 
for even k 
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In Ref. Jy a spin-structure of effective interaction 
within the model of "nearly antiferromagnetic" 
Fermi-liquid was taken into account (spin-fermion 
model). This leads to more complicated combi- 
natorics of diagrams. Spin-conserving part of the 
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Figure 3. Density of states in one - dimensional 
version of the model for incommensurate case and 
different values of T — Vfk/ A. Full lines - ex- 
act numerical diagonalization, dotted lines - our 
Ansatz. 



Figure 4. Energy dependencies of spectral den- 
sity for na = 0.01, t'/t = -0.4, fj,/t = -1.3. 
Incommensurate case. (1) — at the "hot spot" 
p x a/n = 0.142, p y a/ir = 0.857. (2) — close to the 
"hot spot" p x a/n = 0.145, Pyd/n = 0.843. (3)— 
far from the "hot spot" p x a/ir = p y a/iT = 0.375. 



interaction gives formally commensurate combi- 
natorics, while spin- flip scattering is described by 
diagrams with combinatorics of incommensurate 
type. In this case [Q: 

v(k) = ( T f f ° r ° dd \ 0) 
w [ I for even k 

This approach was generalized for calculation 
of two - electron Green's function in Ref. ||, 
where we formulated some recursion relations for 
the vertex part, describing electronic response to 
an external electromagnetic field, allowing calcu- 
lations of conductivity. 

The Ansatz of (|) is in fact not precisely exact 
H , however the analysis of Ref. , shows that it 
is quantitatively good for most interesting cases. 
These conclusions were confirmed in recent pa- 
pers |9|Jl0[|, where the one - dimensional version 
of our model was solved by exact numerical diag- 
onalization. In Fig. 3 we show the comparison of 
results obtained for the density of states in case 
of incommensurate fluctuations by exact numer- 
ics H and via our recursion relations. We can 
see an extremely good correspondence, sufficient 



for any practical purposes. In case of commensu- 
rate fluctuations in one - dimension the Ansatz 
of (Q) misses certain Dyson type singularity of the 
density of states, appearing at the center of the 
pseudogap |) 10|, but this is apparently absent in 
two - dimensional case. 

Consider one-electron spectral density: 
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A(Ep) = ImG R {Ep) 



(10) 



where G R (Ep) is retarded Green's function, ob- 
tained by the usual analytic continuation of (||) 
to the real axis of energy E. In Fig. 4 we show 
typical energy dependencies of A{Ep) obtained 
for the "hot spots" model. More detailed results 
can be found in Refs. Jl|,0|. Similar non Fermi - 
liquid like behavior of the spectral density is eas- 
ily obtained on "hot patches" of the Fermi surface 
shown in Fig.2 (Cf. Ref. g). 

Let us stress that our solution (^) is exact both 
for £ — * cc and £ — > 0, while in the region of finite 
£ it provides apparently very good interpolation. 
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The one-electron density of states: 

N(E)=J2ME,T>) = -lY, ImGR ( E ri (11) 
p p 

is determined by the integral of spectral den- 
sity A(Ep) over the Brillouin zone. Detailed re- 
sults for the "hot spots" model were obtained in 
Ref. |Q , demonstrating smeared pseudogap with 
weak dependence on the correlation length £. For 
"hot patches" model this pseudogap is more pro- 
nounced, depending on the size of these patches. 

"Hot patches" model was applied by us to cal- 
culations of optical conductivity and also (for 
the simplest case of £ — > oo) to the study of the 
main superconducting properties PJl^j. 

Pseudogap phenomena can also be explained 
using the idea of fluctuation Cooper pairing at 
temperatures higher than superconducting tran- 
sition temperature T c . Anticipating the possibil- 
ity of both s-wave and <i-wave pairing, we intro- 
duce the pairing interaction of the simplest (sep- 
arable) form: 



V(p,p') = -Ve(<f>)e(<l>') 



(12) 



where <fi is polar angle determining the direction 
of electronic momentum p in the plane, while for 
e(4>) we assume model dependence: 

f 1 s-wave pairing 

1 \/2cos(2</>) ci-wave pairing 

Analogously to transition from ([!]) to (|^) we in- 
troduce the model interaction (static fluctuation 
propagator of Cooper pairs): 



K//(q)=-AV(£) — 
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(14) 



where A determines the width of superconduct- 
ing pseudogap. We can see that mathematically 
this problem is practically the same as the "hot 
spot" model, but always with combinatorics of 
incommensurate case 0. Then we again obtain 
the recurrence relation for the Green's function of 
the type of (§): 



Sfc(^nCp) 



A 2 e 2 (ct>)v(k) 



ie n - (-l) fc C P + ik(\v x \ + \v y \)n - £ fe+1 (£„£ p ) 



(15) 



where v x = vpcos(p, v y = vpsiiKf), k — £ , 
e n > and v(k) was defined in (||). 

Energy dependencies of the spectral density 
A(Ep) for one-particle Green's function (10), can 
be calculated from (15) for different values of po- 
lar angle cj>, determining the direction of elec- 
tronic momentum in the plane, for the case of 
fluctuations of d-wetve pairing. Calculations show 
that in the vicinity of the point (?r/a, 0) in 
Brillouin zone this spectral density is non Fermi- 
liquid like (pseudogap behavior) . As vector p ro- 
tates in the direction of the zone diagonal the two 
peak structure gradually disappears and spectral 
density transforms to the typical Fermi-liquid like 
with a single peak, which narrows as <f> approaches 
7r/4. Analogous transformation of the spectral 
density takes place as correlation length £ be- 
comes smaller. In case of fluctuation pairing of 
s-wave type the pseudogap appears isotropically 
on the whole Fermi-surface. 
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